Generalized harmonic-fluid approach for the off-diagonal correlations of a 
one-dimensional interacting Bose gas 
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We develop a generalized harmonic-fluid approach, based on a regularization of the effective low- 
energy Luttinger-liquid Hamiltonian, for a one-dimensional Bose gas with repulsive contact interac- 
tions. The method enables us to compute the complete series corresponding to the large-distance, 
off-diagonal behavior of the one-body density matrix for any value of the Luttinger parameter K. 
We compare our results with the exact ones known in the Tonks- Girardeau limit of infinitely large 
interactions (corresponding to K — 1) and, different from the usual harmonic-fluid approach, we 
recover the full structure of the series. The structure is conserved for arbitrary values of the inter- 
action strength, with power laws fixed by the universal parameter K and a sequence of subleading 
corrections. 
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I. INTRODUCTION 

Quasi-one-dimensional quantum fluids display unique 
properties due to the major role played by quantum fluc- 
tuations in reduced dimensionality (see e.g. [![). For ex- 
ample, in an interacting ID Bose fluid quantum fluctu- 
ations destroy the off-diagonal long-range order of the 
one-body density matrix (or first-order coherence func- 
tion), defined as pi(xi,X2) — (ip * (xi)ip(%2)), where ^(a:) 
is the bosonic field annihilation operator at position x. In 
contrast to its 3D counterpart, where the one-body den- 
sity matrix at large distances tends to a constant which 
corresponds to the fraction of Bose-condensed atoms [1] , 
in ID the one-body density matrix at zero temperature 
decays as a power law [E |j| , the coefficient of the power 
law being fixed by the interaction strength: the decay 
is faster as the interaction strength increases from the 
quasi-condensate regime with a fluctuating phase ||, to 
the Tonks- Girardeau regime [6j of impenetrable bosons. 

The one-body density matrix is not only a fundamen- 
tal quantity as it measures the macroscopic coherence 
properties of a quantum fluid with bosonic statistics, but 
also because it is directly related (by Fourier transforma- 
tion) to the momentum distribution of the fluid. While 
the momentum distribution for a Bose fluid is typically 
narrow and peaked around wavevector k = 0, its specific 
form, its broadening due to interactions and its possi- 
ble singularities give a wealth of information about the 
nature of the correlated fluid under study. 

Quasi-one-dimensional Bose fluids find an experimen- 
tal realization in experiments with ultra-cold atomic 
gases confined to the minima of a 2D optical lattice Q. 
For those systems, the momentum distribution is one of 
the most common observables and the one-body density 



matrix has been measured as well [8j , though not yet in 
the quasi- ID geometry. 

From a theoretical point of view, the calculation of the 
one-body density matrix for the ID interacting Bose gas 
has a long history. In the Tonks- Girardeau limit where 
the exact many-body wavefunction is known by means of 
a mapping onto a gas of spinless fermions |(| , the problem 
reduces to the evaluation of a (N — l)-dimensional inte- 
gral. This mathematical challenge has been addressed 
first by Lenard [9fl, and by several subsequent works (see 
e.g. [lO, EE H3.Tl3|)- The main result is the evaluation 
of the large- distance behavior of the one-body density 
matrix in the form of a series expansion, (from [Hj]) 
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where the constant p^ and the coefficients have been 
calculated exactly. In Eq. Jl} and in the following we 
express the one-body density matrix in units of the av- 
erage particle density po and as a function of the scaled 
relative coordinate z = wpo(xi — £2). 

For the case of arbitrary interaction strength, the cal- 
culation of the correlation functions remains a challenge 
(see e.g. [3)) although the model of bosons with con- 
tact repulsive interactions is integrable by the Bcthc- 
Ansatz technique [lq . The power-law decay at large 
distances can be inferred using a harmonic-fluid ap- 
proach |E II EE EE El- The latter is based on an 
effective, low-energy Hamiltonian describing the long- 
wavelength collective excitations of the fluid having a 
linear excitation spectrum u>(k) = v s k. The resulting 
structure for the large-distance series of the one-body 
density matrix reads (from [J]) 
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where K — ■npa/mv s is the universal Luttinger-liquid 
parameter and the coefficients B m are nonuniversal and 
cannot be obtained by the harmonic-fluid approach. By 
noticing that the harmonic-fluid approach is valid also 
in the To nks- Girardeau regime and corresponds to the 
case of Luttinger parameter K — 1, one can directly 
compare the predictions of the two methods. Specifi- 
cally, this comparison shows that the structure of Eq. ([T]) 
is richer than that of Eq. ([2]) obtained by the standard 
harmonic-fluid approach. 

The central result of this work is that by defining a 
properly regularized harmonic-fluid model (to be detailed 
below) the following general structure of the one-body 
density matrix can be obtained for arbitrary values of 
the Luttinger parameter K : 
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where all the coefficients a' n , b m , b' n , c m and c' n are 
nonuniversal and need to be calculated by a fully mi- 
croscopic theory (for possible methods see e.g. [l9l [20|). 
We note that Eq. ^ : (i) generalizes Eq. ^ and (ii) has 
the full structure of the exact result (Q]) in the Tonks- 
Girardeau limit. As a direct consequence of the series 
structure ([3]) the momentum distribution will display sin- 
gularities in its derivatives for k = ±2mirpo. 



II. REGULARIZED HARMONIC-FLUID 
APPROACH FOR BOSONS 

We proceed by outlining the method used. We de- 
scribe the Bose gas with contact interactions (Lieb- 
Liniger model [15[ ) by a harmonic-fluid Hamiltonian, ex- 
pressed in terms of the fields 9{x) and <f>{x) which describe 
the density and phase fluctuations of the fluid [1] : 
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The parameters K and v s entering Eq. (U) above are re- 
lated to the microscopic interaction parameters 2l|, the 
phase field 4>(x) is related to the velocity of the fluid 
v(x) — hV4>{x)/m and the field 6(x) defines the fluc- 
tuations in the density profile. We have adopted here 
an effective low-energy description, which assumes that 
the collective excitations in the fluid are noninteract- 
ing and phonon-like. The description breaks down at 
a length scale a of the order of the typical inter-particle 
distance pQ , as it neglects the broadening and curvature 
of the Bose gas spectrum at finite momentum [TH, [22J. 
Within the harmonic-fluid approximation, the bosonic 
field operator is expressed as ^(x) = \J p(a 



Specifically, its expression in terms of the fields 6{x) 
and 4>{x) reads [3, [2l| 
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where A is a nonuniversal constant, Il(a;) = V9(x)/tt, 
and for compactness of notation we have introduced the 
field 8(x) = 9b + npox + 0(x) 1 where 9b is chosen in 
order to ensure that the average of 9(x) vanish. 

In order to calculate the one-body density matrix, we 
expand the fields 9{x) and <fi(x) in terms of the nor- 
mal modes (bosonic) operators bj, bj which diagonalize 

the Hamiltonian (01 such that Hll = Ylj tiUjbjbj, with 
ujj = v s kj. As we are interested in the thermodynamic 
limit, we have chosen for simplicity the periodic bound- 
ary conditions for a Bose gas in a uniform box of length L, 
where kj = 2nj / L. In this case the mode expansion reads 
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with N and J being the particle number and angu- 
lar momentum operators, and <fio, 9$ being their conju- 
gate fields in the phase- number representation |4]. The 
zero-mode fields </>o, 9q do not enter the calculation of 
the one-body density matrix, as it turns out to depend 
only on the differences 9(x\) — 9{x 2 ) and ^(xi) — 4>{x2). 
In the expressions {6]) and {7} above we have intro- 
duced the short-distance cutoff a ~ p$ , thus regular- 
izing the effective theory. The one-body density matrix 
Pi(xi,X2) — (ip^ {x\)'4>{x2)) is obtained in the generalized 
harmonic-fluid approach from Eq. ([5]) as 

pi(xi,x 2 ) 

= \^\ 2 ([Po + n(xi)] 1/2 e l2m0(xi) e~ l ' t ' {xi) 

(m,m')£Z 2 

x e^^g-^m'e^)^ + n(.x 2 )] 1 / 2 ), (8) 

where the only nonvanishing leading terms satisfy m — 
m'Q. 

We detail now the calculation of the quantum average 
appearing in Eq. ([5]). In order to display its dependence 
only on differences between fields, we re- write the central 
term as 

e i2m©(£Ci) Q -i(j>(xi) e i<f>(x 2 ) g-i2m6(i2) 
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where f(x) = 4>(x) — irx(J)/L and the commutator be- 
tween the 9 and <j> fields is computed in Eq. (fT5j) be- 
low 1 . We then perform a series expansion of the square 
root terms [1 + II/po] 1 ^ 2 in the one-body density ma- 
trix. We define X = H(xi)/p , Y = H(x2)/po, and 
Z = i2m(9(xi) — 9{x2)) ~ i(<p(xi) — ip(x 2 )). Using the 
fact that the fields X, Y and Z are Gaussian with zero 
average, we obtain from Wick's theorem 
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where H n (x) are the Hcrmitc polynomials. To second 
order in X and Y Eq. (flO]) reads 
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N/L = p Q ). Using Eqs. ([T2])-([T4]) we have 
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where a = irpoa. Similarly, the commutator in Eq. @ is 
obtained from Eq. (fTS")) as 
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Notice that the effect of the zero-modes ITo and Jo is ab- 
sent in the thermodynamic limit, because it scales as 1/ L. 
The series expansion in 11/ po is valid for small fluctua- 
tions of the field H(x) compared to the average density p$, 
a/ (X 2 ) < 1 i.e. for a > ^jK/2. By combining the previ- 
ous equations we obtain the final result for the one-body 
density matrix in rescaled units, finding the structure 
displayed in Eq. (J3]): 
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The main expression (|10p requires then the calculation of 
the various two-point correlation functions involving the 
three fields X, Y, Z. All of them can be obtained from 
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(N—(N))/L, and Jo = J— (J). We are now in a position 
to calculate the correlators between the fields X, Y, Z 
in the thermodynamic limit (L — > 00, N — > 00 at fixed 



1 In the following we will not consider the case of a macroscopic 
vorticity in the system, i.e. we assume (J) <C N. 



with Poc = |^| 2 a 1 / 2K c , . 

A few comments are in order at this point. First, 
by taking the limit a — > in Eqs. ([T2|) -(fT5 |) we re- 
cover the results of the standard harmonic-fluid ap- 
proach, i.e. Eq. ([2]). Moreover, the generalized harmonic- 
fluid method produces also the corresponding coeffi- 
cients of the series, namely to the order of approxi- 
mation derived in this work, we have co,o — 1 [1 + 

1 K , K 1. ^ a 2 r , 2K K 2 ]/„ . 

^ (l+«)a 4 r , 4K , K 2 12K+1 K 3 6 l /„ . 

c ,4 - 32Rr 2 [co,o + — + «rrr - ^-if+tJ/ 00 * ' 

Cl)2 ~-2a 2 ^; Cl , 3 ~4a 2 ^+ 1 h- " 

(l+8Jf+4A' 2 )a 2ic+2 ' 



Co,0+2^+if^]/co,o; Ci,4 ~ 

1+4K , K 2 1 I 

t-0,0 -T a i + $k+4K 2 "t" 2a^T+8K+4JC I ' 

1+8K+4K' 2 -1/ c °. 0. However, it should be noted that 
these coefficients do not necessarily coincide e.g. in the 
limit i^T = 1 with the exact ones in Eq. (fTJ) . Our approach 
is still effective, as it suffers from some limitations: first of 
all, we have just used a single-parameter regularization 
which neglects the details of the spectrum of the Bose 
fluid. Second, our approach relies on a hydrodynamic-like 
expression for the field operator ([5]) which implicitly as- 
sumes that the fluctuations of the field U(x) are "small", 
which is not always the case. On the other hand, we see 
from our derivation that the corrections due to the H(x) 
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FIG. 1: One-body density matrix in the Tonks-Girardeau 
limit K — 1 in units of poPoo as a function of the scaled rela- 
tive coordinate z = irpo(xi — X2) (dimensionless) . The result 
of the generalized harmonic-fluid approximation Eq. (|21[l ob- 
tained without taking into account the effect of the field n(:r) 
(solid line) is compared to the exact result Eq. (JTJ (dashed 
line) and to the usual harmonic-fluid approximation Eq. ((2} 
(dotted line), with B = 1, B x = -1/2 and B m>1 = 0. The 
value chosen for the cutoff parameter is a — 1/2. The in- 
set shows the subleading behavior z 1 ^ 2 pi(z) of the one-body 
density matrix in the same notations and units as in the main 
graph. 



fluctuations renormalize the coefficients of the series to all 
orders, giving rise to the contributions in square brackets 
to the Cij above, but do not change the series structure. 

Figure Q] shows our results for the one-body density 
matrix (Eq. (|2ip. solid lines) for K = 1 and a specific 
choice of the cutoff parameter a — 1/2 as suggested by 
the analysis of other correlation functions [23[ . The com- 
parison with the exact result for the Tonks-Girardeau gas 
(Eq. ([1]), dashed lines) yields a reasonable agreement, the 
difference being due to the fact that we have used for the 
sake of illustration the explicit expression for the coef- 
ficients Ci j derived in the current work. It should be 
noted in particular that our generalized harmonic-fluid 
approach restores the correct trend at short distances as 
compared to the usual harmonic- fluid approach (Eq. 
dotted lines). 



III. APPLICATION TO FERMIONS 

It is possible to apply the above approach to the case 
of a ID spinless Fermi gas with odd-wave inter-particle 
interaction [24j, i.e. the ID analog of p-w&ve interac- 
tions. For any value of the coupling strength the (attrac- 
tive) Fermi gas can be mapped onto a (repulsive) Bose 
gas described by the Lieb-Liniger model with dimension- 
less coupling strength given by 73 = — 1/jf [HI, HH, 
the spectrum of collective excitations being the same for 
the two models. Hence, it is possible to describe the 
interacting Fermi gas by the Luttinger-liquid Hamilto- 
nian (H|) as well, except for the so-called Fermi- Tonks- 



Girardeau limit 7f = —00 mapped onto the noninteract- 
ing Bose gas limit where the collective excitation spec- 
trum is quadratic. This allows to estimate the fermionic 
one-body density matrix with a method similar to the 
one described above for the bosons, the fermionic nature 
of the particles being taken into account by a Jordan- 
Wigner transformation Q on the bosonic field opera- 
tor ([5]). In particular, we generalize the series expansion 
known from the usual harmonic-fluid approach, which in 
rescaled units reads (from Q) 
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The coefficients in Eq. (|23j) are nonuniversal and need 
a separate treatment. The limit K = 1 is an excep- 
tion though. It corresponds to the case of noninteracting 
fcrmions where the exact result for the one-body density 
matrix pf^(x 1 ,x 2 ) = J2\k\<k F V , fc( af i)V , *0 c 2), with ipk{x) 
being the single-particle orbitals, in the thermodynamic 
limit and in rescaled units reads 



(0), . sin(z) 
Pif\ z ) = n Po • 



(24) 



Hence, it turns out that for noninteracting fermions all 
the coefficients of the series expression (|2"5|) vanish except 
the leading one. 



IV. CONCLUSION 

In conclusion, we have developed a generalized 
harmonic-fluid approximation and applied it to evalu- 
ate the large-distance behavior of the one-body density 
matrix both for one-dimensional interacting bosons and 
fermions at arbitrary values of the interaction strength 
- the latter expressed in terms of the Luttinger param- 
eter K. In the case of bosons in the Tonks-Girardeau 
limit K=l we recover the full structure of the series ex- 
pansion known for the exact result. In the case of nonin- 
teracting fermions, the exact result shows that the series 
has only one nonvanishing term. In perspective, it could 
be interesting to estimate the coefficients of the series 
expansion at arbitrary interaction strength. 
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